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ABSTRACT 

We describe a coagulation model that leads to the rapid formation of super- 
Earths and the cores of gas giant planets. Interaction of collision fragments 
with the gaseous disk is the crucial element of this model. The gas entrains 
small collision fragments, which rapidly settle to the disk midplane. Protoplanets 
accrete the fragments and grow to masses > 1 in ~ 1 Myr. Our model 
explains the mass distribution of planets in the Solar System and predicts that 
super-Earths form more frequently than gas giants in low mass disks. 

Subject headings: planetary systems - solar system: formation - planets and 
satellites: formation 



1. INTRODUCTION 



Collisional cascades play a central role in planet formation. In current theory, planets 
grow from collisions and mergers of km-sized planetesimals in a gaseous disk. As planets 
grow, they stir leftover planetesimals along their orbits to high velocities. Eventually, col- 
lisions among planetesimals produce smaller fragments instead of larger, merged objects. 
Continued stirring leads to a cascade of destructive collisions w hich grinds the leftovers to 
dust. This process (i) ex plains the masses of terrestrial planets (IKenyon fc Bromley II2006I ) 
and Kuiper belt objects (IKenyon et al. 1120081 ) and (ii) prod uces debris disks similar to those 
observed around nearby main sequence stars (IWyatt II2008I ). 
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Numerical simulations of icy pl anet formation suggest the cascade limits the masses of 
growing protoplanets to ~ 0.01 Mq (IKenyon fc Bromley II2008I . hereafter KB08). This mass 
is much smaller than the core mass, > 0.1-1 M^, r equired for a protoplanet to accrete ga s 



and become a gas giant planet fjPoUack et al. Ill996l : llnaba et al. II2003I : lAlibert et al. II2005I ). 



Unless icy protoplanets can accrete collision fragments before the fragments are ground to 
dust, these protoplanets cannot grow into gas giant planet cores. Thus, finding a mechanism 
to halt the cascade is essential to form gas giant planets. 

Here, we describe how interactions between the fragments and the gaseous disk can halt 
the cascade. In our picture, the gas traps small fragments with sizes of 0.1 mm to 1 m and 
prevents them from colliding at large velocities. These fragments then settle rapidly to the 
disk midplane, where protoplanets can accrete them. For a broad range of initial conditions, 
analytic results and detailed numerical simulations demonstrate that this process yields 1- 
10 Me cores in 1-2 Myr. 

We develop the analytic theory in §2 and derive the conditions needed for protoplanets 
to accrete collision fragments and grow to masses of ~ 1 in 1-2 Myr. We confirm these 
estimates in §3 with detailed numerical calculations. We conclude with a brief discussion in 
^4. 



2. PHYSICAL MODEL 



The crucial element of our model is the interaction of collision fragments with the gaseous 
disk. Fragments larg er than the 'stopping radius' ^ 0.5-2 m at 5-10 AU (jWeidenschilling 



19771 : iRafikov I I2004J ) . orbit with the growing protoplanets, independently of the gas. De- 
structive collisions among these fragments fuel the coUisional cascade. However, the gas 
entrains particles with radii r < r^. These fragments orbit with the gas; thus, their velocity 
dispersions are small and independent of massive protoplanets. By trapping small collision 
fragments, the gas halts the coUisional cascade. 

The gas also allows protoplanets to accrete the debris. When the coUisional cascade 
begins, the mass in leftover planetesimals is ~ 1-10 M®. The cascade grinds all of this 
mass into small fragments which are trapped by the gas. Most of the trapped fragments fall 
through the gas into the midplane of the disk, where growing protoplanets a ccrete them. Pro 
toplanets that accrete ~ O-ItI before the gas dissipates ( ~ 3-10 Myr; iHartmann et al. 



19981 : iHaisch. Lada. &: Lada Il200ll : Kennedy fc Kenyon II2009I ) become gas giants. Thus, our 
model yields gas giant cores if (i) the coUisional cascade produces fragments fast enough, 
(ii) the fragments quickly settle to the midplane, and (iii) the largest protoplanets rapidly 
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accrete the fragments. 



To examine whether this physical model leads to cores with masses of ~ 1 M^, we 
consider the growth of planets in a disk of gas and icy objects around a star with mass 
M^. Material at a distance a from the central star orbits with angular frequency VL and 
has surface densities (solids) and Eg (gas). We adopt a solid-to-gas ratio of 1:100 and 
Tig = T,sfl Xm ot'^^'^ 1 where S^^o = 2.5 g cm~^ at 5 AU and Xm is a scale factor. 

Forming icy protoplanets is the first step in our model. In an ensemble of 1 km planetes- 
imals, collisional growth yields a few 100 km objects - ' oligarchs' - that contain an ever- 



increasiiig frac tion of the mass in solids (llda &: Makino I Il993l : IWetherill &: Stewart I Il993 



Rafikov II2003I ). From numerical simulations of planet growth at 30-150 AU, the timescale 
to produce an oligarch around a solar-type star is (KB08) 



^1000 



10^ xZ}-^"" 



5 AU 
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(1) 



Thus, oligarchs form at 5 AU before the gas dissipates. 



Once oligarchs form, collisions among leftover planetesimals produce copious amounts 
of fragments. In the high velocity limit, the collision time for a planetesimal of mass M in a 
swarm of icy planetesimals with ma ss M, radius r, and surface density S is tc ~ (S vr fi) 
JColdreich. Lithwick. fc Sari]l2004h . Thus, 
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(2) 



Collisions among planetesimals produce debris at a rate M ^ N 6M t~^, where N is the 
number of planetesimals of mass M and 6M is the mass in fragments produced in a single 
collision. In an annulus of width 6a at distance a from the central star, N ^ 27t a 6a/ M. 
If ~ 10% of the mass in each collision is converted into fragments 
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(3) 



where we have set t he width of the annulus equal to the width of the 'feeding zone' for a 
0.1 protoplanet (ILissauer 1119871 ). Thus, disks with Xm ^ 1-2 produce fragments at a 
rate sufficient to form > 1 cores in 1-2 Myr. 

Most of the mass in fragments settles quickly to the disk midplane. For a settling time 
ts ~ 4a;m(l m/r) yr (IChiang fc Goldreich 1 119971 ). fragments with r > 0.1 mm reach the 



midplane on the colli sion timescale of ~ 10^ yr (Eg. [21). For a siz e distribution n(m) cx m'^ 
with g = -1 to -0.8 JPohnanvi II1969I : Iflolsapple fc Housen |[2007h . 66% to 99% of the total 
mass in fragments with r < 1 m settles to the midplane in < 10^ yr at 5-10 AU. 
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Oligarchs rapidly accrete fragments in the midplane. The maxi mum accretion rate for 



an ohgarch with ~ 0.01 at 5 AU is ~ 5 x lO'^ yr"^ flRafikov I l2004f ). This 



maximum rate yields 5 cores in 1 Myr. At the onset of the cascade, our simulations 
suggests oligarchs at 5 AU accrete at a rate 



M„ ^ 10 



-6 



2/3 



M 



f 



Ma 



yr 



^0.01 M®; VI 
where Mf is the total mass in fragments in a feeding zone with width 6a 
Thus, protoplanets likely reach masses ~ 1 in 1-2 Myr. 



(4) 

0.2 AU at 5 AU. 



These analytic estimates confirm the basic aspects of our model. In a gaseous disk with 
250 g cm~^, the gas halts the coUisional cascade. Collision fragments entrained by the 
gas rapidly settle to the midplane. Protoplanets with masses ~ 0.01 can accrete collision 
fragments rapidly and grow to masses ~ 1 Mq before the gas dissipates in 3-10 Myr. 



NUMERICAL MODEL 



To explore this picture in more deta il, we calculate the formati on of cores with our hybrid 
multiannulus coagulation-n-body code (IBromley fc Kenyon II2006I ). In previous calculations, 
we followed the evolution of objects with r > r^; collision fragments with r < were 
removed by the collisional cascade (KB08). Here, A ve inclu de the evolution of small fragments 
entrained by the gas. We follow iBrauer et al. I (j2008al ) and calculate the scale height of 
small particles wit h r < r.^ as if = aHgl\Tc im.{St,Q.b){l + St)], where Hg is the scale 
height of the gas (IKenyon fc Hartmann 1 119871 . KB08), a is the turbulent viscosity of the 
gas, and St = rpsQ/csPg. In this expression for the Stokes number (St), Cg is the sound 
speed of the gas, pg is the gas density, and ps is the mass density of a fragment. We 
assume small particles have vertical velocity v = HQ and horizontal velocity h = l.Qv. 
Protoplanets accrete fragments at a rate naVrei, where n is the number density of fragments, 
a is the cross-section (including gravitational focusing), and Vrei is the relative velocity (e.g., 
Kenyon fc Luu 1119981 . Appendix A. 2). Although this app roximation neglects many details of 
the motion of particles in the gas (IBrauer et al. Il2008al ). it approximates the dynamics and 
structure of the fragments reasonably well and allows us to calculate accretion of fragments 
by much larger oligarchs. 



Using the statistical approach of ISafronov I (119691 ). we evolve the masses and orbits of 
planetesimals in a set of concentric annuli with widths Sai at distances from the central 
star (KB08). The calcul ations use realistic cross-sections (including gravitational focusing) 
to derive collision rates ( Spaute et al. 1 1991) and a Fok ker-Planck algorithm to derive grav- 



itational stirring rates (lOhtsuki. Stewart, fc Idall2002l ). When large objects reach a mass 
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Mproi we 'promote' them into an n-body code (IBromley &: Kenyon II2006I ). This code follows 



the trajectories of individual objects and includes algorithms to allow interactions between 
the massive n-bodies and less massive objects in the coagulation code. 

To assign collision outcomes, we use the ratio of the center of mass collision energy Qc 
and the energy needed to eject half the mass of a pair of collidin g planetesimals to infinity 
Q*^. We adopt Q*^ = Qbrf^" + QgpgV^'' Jsenz fc Asphaug I fl999l ). where Qbr'^" is the bulk 
component of the binding energy, QgPgV^o is the gravity component of the binding energy, 
and pg is the mass density of a planetesimal. The mass of a merged pair is M]_+JW2_^^il^ 
wher e the mass ejected in the collision is Mej 
1999h . 



0.5(Mi + M2)iQc/Q*dr^ flKenvon fc Luu 



Consistent with recent ra-body simulations, we consider two sets of fragmentation pa- 
rameters fi. Strong planetesimals have = (Qb = li lO^ i or 1 0^ erg g~^, ~ 0, Qg 



= 1.5 erg g' 
have fw = {Qb 



1 1.25 



Leinhardt fc Stewart 



2 X 10^ 
2008h . 



1.25; KB08, iBenz fc Asphaug Ill999l ). Weaker planetesimals 



— 1 4 

erg g cm , 



-0.4, Qg = 0.22 erg g-^ cm'^-^ I3g = 1.3; 



Our initial conditions are appropriate for a disk aroun d a young star (e.g. lDuUemond fc Dominik 



2005; 


Ciesla 


2007a: 


Garaud 


2007; 



with Oj = 5-10 AU and Stti/ai 
radii ranging from r^m = '"s " 
and an initial mass distributio n riiiMik) oc Mt ' 
5 = Mik+i/Mik 
eccentricity eo = 



= 0.025. The disk is com posed of small planetesimals with 
0.5-2 m jRafikovl[20oi ) to tq = 1 km, 10 km, or 100 km 



0.17 



[e.g., 



Kenyon fc Luu 



= 1.4-2 
10~^ and inclination io 



The mass ratio between adjacent bins is 
1998. KB08). Each bin has the same initial 



eo/2. 



For each combination of tq, fi, and Xm = 1^5, we calculate the growth of oligarchs with 
two different approaches to grain accretion. In models without grain accretion, fragments 
with r < rmin are 'lost' to the grid. Oligarchs cannot accrete these fragments; their masses 
stall at M < 0.1 M®. In models with grain accretion, we track the abundances of fragments 
with 0.1 mm < r < rmin which settle to the disk midplane on short timescales. Oligarchs 
can accrete these fragments; they grow rapidly at rates set by the production of collision 
fragments. 



For the gaseous disk, we adopt a = 10 an initial surface density, Eg^o = 100 Ss,o x 



and a depletion time tg = 3 Myr. The surface density at later times is S„ 



'9,0 



-t/tg 



(Lin & Papaloizou 


1986; 


Ward 


1997 


)• 



growth of protoplanets; thus our approach underestimates the growth time. We also ignore 
the radial drift of fragments coupled to the gas. Depending on the internal structure of 
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the disk, fragments can drift inward, drift outward, or become concentrated within local 



pressure maxima or turbulent eddies ( 


Weidenschi 


line 


1977 




Haehiehipour & Boss 


2003; 


Inaba & Baree 


2006: 


Masset et al. 


2006 




Ciesla 


20071 


); 


Kretke & Lin 


2007 




Kato et al. 



20081 ). Here, our goal is to provide a reasonable first estimate for the growth rates of proto- 



planets. We plan to explore the consequences of radial drift in subsequent papers. 



4. RESULTS 

Fig. [1] shows mass histograms at 1-10 Myr for coagulation calculations without grain 
accretion using vq = 1 km and the strong fragmentation parameters [fs)- After the first 
oligarchs with M ~ 0.01 form at ~ 0.1 Myr, the collisional cascade starts to remove 
leftover planetesimals from the grid. Independent of Qb, the cascade removes ^ 50% of the 
initial mass of the grid in ~ 4 Myr. As the cascade proceeds, growth of the largest 

oligarchs stalls at a maximum mass Mo^max ~ 0.1 Xm Mq. 

These results depend weakly on tq. The time to produce the first oligarch with r ~ 
1000 km increases with tq, tiooo ~ 0.1 x^^'^^ {ro/10 km)^/^ Myr. Calculations with larger tq 
tend to produce larger oligarchs at 10 Myr: Mo^max ~ 1 M® (2 M®) for tq = 10 km (100 km). 
In ~ 50 calculations, none produce cores with Mo^max ^ 1 M^ on timescales of < 10 Myr. 

For ro < 100 km, our results depend on /j. In models with ro = 1 km and 10 km, the 
fw fragmentation parameters yield oligarchs with smaller maximum masses, Mo^max ~ 0.3- 
0.6 Mq. Because leftover planetesimals with r ~ 1-10 km fragment more easily, the cascade 
begins (and growth stalls) at smaller collision velocities when oligarchs are less massive 



(IKenvon et al. Il2008[ ). 



Calculations with grain accretion produce cores rapidly. Fig. [2] shows results at 1- 
10 Myr for calculations with tq = 1 km and the fs fragmentation parameters. As the first 
oligarchs reach masses of ~ 0.01 M^ at 0.1 Myr, the cascade generates many small collision 
fragments with r ~ 1 mm to 1 m. These fragments rapidly settle to the disk midplane and 
grow to sizes of 0.1-1 m. When the cascade has shattered ~ 25% of the leftover planetesimals, 
oligarchs begin a second phase of runaway growth by rapidly accreting small particles in the 
midplane. For calculations with Xm = 1^5, it takes ~ 1-2 x^''^^ Myr to produce at least one 
core with Mo > 1-5 Mq. Thus, cores form before the gas dissipates. 

These results depend on tq. For ro = 10 km, fragmentation produces small grains 2-3 
times more slowly than calculations with ro = 1 km. These models form cores more slowly, 
in 5-10 Myr instead of 1-2 x^'"^^ Myr. For models with ro = 100 km, fragmentation 

yields a negligible mass in small grains. Thus, cores never form in < 10-20 Myr. 
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The timescales to form cores also depend on /j. Calculations with the fw parameters 
form cores 10% to 20% faster than models with the fs parameters. 



CONCLUSIONS 



Gaseous disks are a crucial element in the formation of the cores of gas giant planets. 
The gas traps small collision fragments and halts the coUisional cascade. Once fragments 
settle to the disk midplane, oligarchs accrete the fragments and grow to masses > 1 in 
1-3 Myr. 

Our model predicts two outcomes for icy planet formation. Oligarchs that form before 
(after) the gas disk dissipates reach maximum masses > 1 M® (< 0.01-0.1 M®). Setting the 
timescale to form a 1000 km oligarch (Eq. [1]) equal to the gas dissipation timescale tg yields 
a boundary between these two types of icy protoplanets at ^ 15 x^'^(tg/3 Myr)^/^ AU. 
We expect massive cores at a < and low mass icy protoplanets at a > a^. 

This prediction has a clear application to the Solar System. Recent dynamical calcula- 
tions suggest that the Solar System formed with fou r gas giants at 5-15 AU and an ensemble 
of Pluto-mass and smaller objects beyond 20 AU (IMorbidelli et al. II2008I ). For a protoso- 



lar disk with xl^ {tg/'i Myr) ^ 1, our model explains this co nfiguration. Disks with thes e 



parameters are also common in nearby star-forming regions (lAndrews fc Williams II2005I ). 
Thus, our results imply planetary systems like our own are common. 

Our model yields a large range of final masses for massive icy cores. Protoplanets that 
grow to a few well befor e the gas dissipates can accrete large amo unts of gas from the 



disk and become gas giants (jPoUack et al. Ill996l : lAlibert et al. II2005I ). Protoplanets that 



grow m ore slowly cannot accrete much gas and become ic y 'super-Earths' with much lower 



masses ( iKennedy et al. I l2006l : Kennedy fc Kenyon II2008I ). For solar-type stars with t 



3 Myr, our results suggest that gas giants (super-Earths) are more likely in disks with Xm ^ 
1-5 (xra < 1.5) at 5-10 AU. 

Testing this prediction requires (i) extending our theory to a range of stellar masses 
and (ii) more detections of massive planets around lower mass stars. We plan to explore 
the consequences of our model for other stellar masses in future papers. Larger samples of 
planetary systems will test the apparent trend that g as giants (super-Earths ) are much more 



common around solar-type (lower mass) stars (e.g.. ICumming et al. II2008I : iForveille et al 
20081 ). Comparing the results of our planned numerical calculations with these additional 
observations will yield a clear test of our model. 
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Fig. 1. — Mass histograms at 1 Myr (left panel), 3 Myr (center panel), and 10 Myr (right 
panel) for coagulation calculations without grain accretion using the fs fragmentation pa- 
rameters at 5 AU. Magenta histograms plot median results for 25 calculations with Xm = 1; 
blue histograms show median results for 25 calculations with Xm = 5. Independent of disk 
mass, calculations without grain accretion yield planets with maximum masses < 1 M^ in 
10 Myr. 
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Fig. 2. — As in Fig. 1 for calculations with grain accretion. When large oligarchs can accrete 
fragments trapped by the gas, disks with Xm > 1 produce gas giant cores in 3-10 Myr. 



